Collisions, Cosmic Radiation and the Colors of the Trojan 

Asteroids 

M.D. Melita^ 

Institute de Astronomia y Fisica del Espacio. Buenos Aires. Argentina. 

me 1 i t aO i af e . ub a . ar 

and 

G. StrazzuUa^ 
INAF-Osservatorio Astrofisico di Catania, Italy 

gianniOoact . inaf . it 

and 
A. Bar-Nun^ 

Department of Geophysics and Planetary Sciences. Tel Aviv University, Tel Aviv,Israel. 

aJsivabSpost . tau .ac.il 

Received ; accepted 



Submitted to ICARUS 



- 2 - 



ABSTRACT 



The Trojan asteroids orbit about the Lagrangian points of Jupiter and the 
residence times about their present location are very long for most of them. 
If these bodies originated in the outer Solar System, they should be mainly 
composed of water ice, but, in contrast with comets, all the volatiles close to the 
surface would have been lost long ago. Irrespective of the rotation period, and 
hence the surface temperature and ice sublimation rate, a dust layer exists always 
on the surface. We show that the timescale for resurfacing the entire surface of 
the Trojan asteroids is similar to that of the flattening of the red spectrum of 
the new dust by solar-proton irradiation. This, if the cut-off radius of the size 
distribution of the impacting objects is between 1mm and Im and its slope is -3, 
for the entire size-range. Therefore, the surfaces of most Trojan asteroids should 
be composed mainly of unirradiated dust. 

Subject headings: Trojan Asteroids. Asteroids, Surfaces. Comets. Solar System, 
Origin. 
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1. Introduction 

All the Trojan asteroids orbit about the Sun at roughly the same heliocentric distance 
as Jupiter, residing close to the Lagrangian stable points, leading or trailing at ~ 60° from 
the planet. They are particularly interesting because of being relatively isolated, lying just 
beyond the snow line, the distance at which water can exist in the form of ice. As such, it 
may be expected that the bulk of the Trojans were mainly composed of ices. However, at 
the moment, no water or other more volatile substances have been detected on the surface 
of a Trojan (Yang & Jewitt 2007, Emery & Brown 2003; 2004). In the case of Ennomos, an 
unusually hight-albedo object, the surface-content of water ice has been quantified to be 
below 10% in mass (Yang & Jewitt 2007). However the lack of detection does not imply 
that water ice is not there. Brunetto and Roush (2008) have in fact shown that a relatively 
thin (tens of microns) crust of specific refractory materials (they considered irradiated 
methane ice) can mask the presence of water ice bands in the spectrum. 

The Trojans have been characterized as similar to cometary nuclei or extinct comets 
(Jewitt and Luu 1990). Most of the objects fall into the D and P asteroid taxonomic classes 
and a few objects show spectra with neutral or even negative slopes in the visual range 
(Bendoya et al. 2004), with the majority having moderately-red surfaces (see figured]). 
They have been observed to posses low albedos (Fernandez et al. 2003) of about 4% and 
infrared observations indicate the presence of silicates (Cruikshank et al. 2001, Emery et al. 
2006) and depletion of water ice and volatiles (Emery & Brown 2003; 2004). 

Numerical models of the collisional evolution of the Trojan asteroids (Marzari et al. 
1997) are successful in explaining the present Trojan population (Jewitt et al. 2000). The 
steep slope of the Trojan size distribution at diameters larger than about 50 — 100km is 
essentially unaltered during the collisional evolution and must reflect the formation process 
of these bodies. At smaller sizes, collisions have produced a power law size distribution. 
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with a slope characteristic of collisionally relaxed populations. The formation of prominent 
dynamical families in the Trojan swarms is a natural outcome of the collisional process, and 
indeed a significant number of families are known (see for example Beauge & Roig 2001). 

It has also been found that the absolute-magnitude distribution of the Trojans in 
stable orbits is bimodal, while that of the unstable orbits is unimodal, with a slope similar 
to that of the small stable Trojans (Melita et al. 2008a). This supports the hypothesis that 
the unstable objects are mainly byproducts of physical collisions, the transitional size has 
been found to be ~ 9.5. But the slopes of the reflectance spectra show no relation with 
the size of the objects (see figure or with the stability of the orbits (Melita et al. 2008a). 

We intend to understand why most of the Trojans are moderately red, regardless if 
they are large bodies or fragments of collisions. These asteroids are subject to impacts 
which can extract fresh material from the interior and scatter it on the surface. The fresh 
material would be composed of dust and ice. After the Deep- Impact event, we note the 
increase of emission of both organic compounds and water ice (Lisse et al. 2006). But, if 
water ice were deposited on the surface after an impact, it would sublimate very rapidly. 
Therefore, after a physical impact, a layer of unirradiated ice -which has a typically red 
spectral slope- is deposited in the surface. It has been suggested that terrestrial bitumens 
are good spectral analogues for such a dust layer (e.g. Moroz et al 2004 and references 
therein). The less energetic solar protons modifies the spectroscopic properties of this 
dust, reducing its albedo and flattening the spectral slope (Moroz et al. 2004). But if the 
collisional process can resurface the asteroids on a timescale which is similar, or shorter, 
than the timescale that flattens the spectrum, it is expected that most of the Trojans will 
be observed to be red. In the following section we present a model to estimate relevant 
timescales and we discuss its results in the last section. 
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2. Surface alteration timescales 

2.1. Timescale of flattening of the spectral slope by irradiation 

There are a number of energetic ion populations bombarding the surfaces of small 
bodies in the outer Solar System: solar wind ions, ions from solar flares, galactic cosmic 
rays, and the so called "anomalous" cosmic ray component. The last component is however 
relevant only to objects at the Pluto orbit and beyond. 

Solar wind is an expanding flux of fully ionized plasma that reaches, at distances 
greater than a few solar radii, an expansion speed of about 400 km s"^ corresponding to an 
energy of R^l keV/amu (slow solar wind, see e.g. Goshng 2007). At 1 AU its flux is r^IO^ 
protons cm'^s^^ and it decreases with the square of the solar distance. The effects of the 
solar wind ions irradiating "red" organic bitumens (namely asphaltite and kerite) have been 
simulated in the laboratory (Moroz et al 2004). The hypothesis is that terrestrial bitumens 
are good spectral analogues of the red organic surfaces of small objects in the Solar System 
(we consider Trojans in the present paper). The main result is a progressive flattening 
of the red spectrum of bitumens with ion fluence. Moreover Moroz et al (2004) suggest 
that the mechanism that produces the flattening is connected to the energy deposited by 
elastic collision of the incoming ions with the atoms in the target. This allows to extend the 
laboratory results obtained with lO's keV H+, N+, and Ar"*""*" to all solar wind ions. Because 
in the laboratory the flattening of the asphaltite spectrum in the 0.3-0.8 micrometers 
spectral region has been obtained after about 0.4 x 10^^ C-displacements cm~^ (see Moroz 
et al, top panel in Fig 14), we need to calculate the time necessary to solar wind ions (at 
5 AU) to cause the accumulation of such an amount of damage on the bombarded surface. 
As said, solar wind ions have energy of about 1 keV/amu. Most of them are protons (96 %) 
and He (4 %) with traces of heavier elements. Using the TRIM code (Transport of Ions 
in Matter; e.g., Ziegler et al. 1996) we calculate that the number of displaced atoms in 
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asphaltite is about 1.5 for 1 keV protons and 20 for 4 keV alpha particles. Considering the 
respective abundances, we obtain an average of 2.2 C-displacements per solar wind ion (the 
contribution of heavier ions can be neglected) . At a distance of 5 AU from the Sun we then 
expect that an exposed surface will suffer about 10^ x 2.2/25= 9x 10^ C-displacements 
cm~^sec~^. Thus the time necessary to reproduce the effects observed in the laboratory is: 
0.4 X 10^79 xl06= 5.5 X 10^° sec i.e. about 2x10^ yrs. 

It is important to note that solar wind particle irradiation affects only a thin surface 
layer (< 2000 A). This is however sufficient to affect the optical properties of the surface 
(as simulated in the laboratory), but may be easily destroyed by impact resurfacing. As 
discussed here, such resurfacing of irradiated material probably occurs constantly. 

In order to produce a radiation-induced damage of deeper layers, more energetic ions 
should be considered. Solar cosmic rays are produced during solar flares. An average flux 
of solar protons having an energy of 100 keV and a penetration depth of about 1 fim has 
been estimated to be a; 2x10^^ protons cm~^ yrs~^ (Lanzerotti et al., 1978). It is inversely 
proportional both to the square of the energy and to the square of the solar distance. At 5 
AU the time necessary to reach the maximum C-displacements of our experiments is about 
2x 10^ yrs. 

Galactic cosmic rays are less abundant but significantly more penetrating. In a recent 
review Hudson et al. (2008, and references therein) reported the estimate of the radiation 
doses for objects at different distances from the Sun (see their Table 1). It has been 
found that objects orbiting at 5-35 AU accumulate in about 10^ yrs an irradiation damage 
sufficient to be fully chemically altered down to depth of the order of 1 meter. Although 
galactic energetic particles transfer most of their energy through electronic loss (ionizations 
and excitations), however at the end of their path they slow down and are thought to 
produce, in the so called Bragg peak, the same effects meassured in the laboratory for low 
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energy particles. 

In summary we have shown that the uppermost layers of an object orbiting at the 
Jupiter distance can be weathered in about 10^ yrs. If its surface is an organic material 
spectrally similar to bitumens, its red sloped spectrum is flattened on that time scale. 
Deeper layers are weathered in longer times: at a depths of the order of 1 meter the target 
is fully altered in about 10^ yrs. 



2.2. Timescale of sublimation of amorphous water-ice 

The timescale for subhmation-mantle growth, tm, is (Jewitt 2002): 

P L 
Jm dm/dt^ 

where L is the depth of the sublimated material, p is the mass density of the asteroid, /m 
is the fraction of the sohd mass that cannot be ejected by gas drag and dm/dt is the drag 
mass-loss rate per unit area. 

For a fast rotator, the effective temperature at the surface of a Trojan is ~ 120°K 
-assuming an albedo of 0.04. At this temperature, the mass loss rate of perfectly absorbing 
water-ice in equilibrium sublimation is approximately dm/dt w 3 lQ~^^kg s~^. To find 
a reliable estimation for dm/dt at T — 120° K is quite a complicated problem and we have 
reached this value after much deliberation, following the comments of the reviewer. The 
water-vapor pressure in torr, P, vs. temperature, T, was calculated, using the polynomial 
from the International Critical Tables between 140"^ and 1S0°K: 

log (P) = -2445^-5646 ^ g 23;l21 x log(r) - 0.167706 xT+ 1.20514 x lO'^T^ - 6.757169, 

and checked against the P vs. T in the CRC-Handbook of Chemistry and Physics and 
found to agree with Sack and Baragiola (1993). The pressure at 120°K is 1.8 x 10~^° Pa. 
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Since the flux of sublimating water molecules, F, is, 



(27rXmT)V2' 

where K is Boltzmann's constant and m the mass of the molecule. We recalculated the square 

root term in the denominator and found it to be: {27rKmTy^'^ — 1.77 x lO'^^kg 

per molecule. Hence the flux calculated from the pressure as obtained by the polynomial is: 

PiT = 120"K] 

F(T = 120°X) = , = 1.1 X lO^^molec. s'^ = 1.1 x 10^ molec. cm~^ 

This flux was compared with the flux meassured by Ntesco and Bar-Nun (2005) in their 
Figure 1 and we reahzed that the measured value at 120° K is the background of water 
vapor in our vacuum chamber. Extrapolating the flux vs T plot in their Figure 1 to 120"i^' 
leads to a flux of approximately 10^^ molec. cm~^ s~^, which is 100 times larger than the 
one calculated by the P vs T through the polynomial, due to the very large surface area 
of the amorphous ice, much larger than the area calculated by the geometric dimension of 
the substrate on which the ice was deposited. As measured by Bar- Nun et al. (1987), by 
Ar adsorption at 120°K, a surface area of 38m^ was found, in good agreement with 
Mayer and Pletzer (1986) who measured 40m^ g~'^ at 11^°K. Sack and Baragiola' s (1993) 
Figure 2 data can be extrapolated to 8 x IQ^^molec. cm^^ at 120"A' and shows an 8 
fold increase in the vapor pressure in amorphous non annealed ice, which they attribute to 
the larger surface area. These authors (Figure 3)show a 7 fold increase in the flux at 135"^ 
with an ice thickness between 0.12 and 5.1 /xm and a 7 — 8 fold decrease at 135"^ during 
200mm at that temperature. Both observations point to a considerable effect in surface 
roughness. To be on the conservative side we shall adopt a flux of 5 10^° molec. cm~^ or 
1.5 10~^^kg m"'^ s"^. It should be remembered that slight deviations from an albedo of 0.4 
or the fast rotator approximation could result in large changes in the rate of sublimation. 
The question is then which flux of water molecules from an ice surface at 120° K should 
be adopted?. The one obtained from Notesco and Bar-Nun's (2005), 10^^ molec. cm""^ 
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or Sack and Baragiola's (1993) extrapolated 8 x 10^° molec. cm~'^ s~^7. Since this mass 
loss rate is very low, even the smallest sub- micron sized grains on the surface of the Trojan 
asteroids are retained, hence, for the fast rotator case, we take Jm = 1- 

The known rotation periods of Trojans range approximately from 4hr to 40/ir 



(Hartmann et al. 1988, Binzel and Sauter 1992, Melita et al. 2008b 

As one extreme case, in the fast rotator approximation, the time that it takes to 
sublimate a layer equal to the penetration-depth of the solar protons, roughly a lOOyum dust 
layer, is only of ~ 200yr. 

In the case of the slower rotators we may take, as another extreme case, the stationary 
approximation. The stationary Trojan temperature of the side facing the Sun can be 
estimated using the equation of effective temperature with back irradiation from half a 
sphere, rather than the entire surface of the sphere as in the case of a fast rotator. This 
leads to a higher temperature, namely Ts ~ 1Q8°K, and the corresponding mass-loss rate 
dm/dtiTs) ^ 3 IQ-^kg mr'^ s^^ (CRC-Handbook of Chemistry and Physics 1981-1982). In 
this case, the fraction of the solid mass that cannot be ejected by gas drag can be estimated 
as (Jewitt 2002): 

_ log^au/cic) 

JM 



log{aM/am) 

where and the maximum and the minimum grain sizes in the distribution 

respectively and Oc, the maximum particle-radius that is retained, given by: 

Ve{R) dm/dt{Ts) 

""'^ 'gT? ' 

where Ve{R) is the escape velocity from a body of size R and G is the gravitational constant. 
If we take Oa/ = 1mm and a.^ = Ifim and R = 100km, then Oc ~ 2/xm. We find that, for a 



^Rotation periods of Jupiter Trojan asteroids, determined before March 2006, can be 



found at the Minor Planet Center website, http://cfa-www.harvard.edu/cfa/ps/mpc.html 
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body of rotational period P — 40/ir, the timescale in which an insulator-mantle of thickness 
L — lOO/im grows, is only tm{P — 40/ir) ~ Sdays. 



2.3. Timescale of collisional re-surfacing 

To compute the timescale of the collisional resurfacing process we have adapted the 
model by Gil-Hutton (2002) and incorporated the latest knowledge gained by the Deep 
Impact mission (see for example Richardson et al. 2007). 

The time to cover the entire surface of an asteroid of radius R with debris extracted 
from below the surface by the action of physical collisions, tcr, is given by: 

1 

Of 

where Sr — S/{At^R'^) is the fraction of the surface that is covered by material excavated by 
collisions per year and S is the total area covered by coUisional ejecta each year. 

The total area covered by ejecta per year, S, for an object of radius R is: 

N{t) Ae{t) dr, (1) 

n 

where N{r) is the number of collisions per unit time and AE{r) is the area covered by the 
ejecta produced when projectiles of radii r collide with an object of radius i?, at a typical 
encounter velocity v. We write Ae as Ae = I x Ac{R,r,v), where Ac{R,r,v) is the area of 
the crater produced and A^ = t: d"^, where d is the diameter of an idealized crater, and we 
assume the densities of the target and the impactor to be equal, i.e. PR — Pr — l-Bgcm'^. 
The volume of the crater in the gravity dominated cratering regime is estimated as 
(Richardson et al. 2007, Holsapple 1993): 
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where g{R) is the surface gravity of the target, m is the mass of the impactor, v its impact 
velocity, p is the density which is assumed to be equal for the target and the impactor. The 
scaling constants Ki = 0.24 and fi = 41, correspond to a loosely bound material such as 
sand, which gave a good agreement with the size of the crater produced in the Deep impact 
mission (Richardson et al. 2007). The diameter of the crater is easily obtained from 

Now, we can write N as, 

N = Pi {R + rf dN{r), 

where Pi is the intrinsic probability of collisions, as defined by Wetherill (1967), is the 
probability of collision per unit time, per unit cross section and per number of colliding 
pairs and dN{r) is the number of objects with radii r. The values of Pi and v, have 
been computed specifically for the Trojans by Dell'Oro et al. (1998) giving approximately 
6.5 10^"*^^ km~'^ yrs~^ and 5 km respectively. 

The integral on equation [1] depends critically on the limits, which must be handled 
carefully. The upper limit, Tmax-, is the maximum radii of the projectiles that do not shatter 
a target of radius R. For the p/j = pr case is, 

W = (5 ^vl - l)-l/^ 

(Petit & Farinella 1993), where (3 = 10~^ g erg~^ is a crater excavation coefficient. While 
rmin is the cut-off radius. 

Finally, to compute tc, we need an estimation of the observed distribution of sizes in 
the Trojan population, which, for the smallest objects is (Jewitt et al. 2000), 

dN{r) = 1.5 10' ("i^y. (3) 
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2.3.1. The choice ofrmin, the cut-off radius 

The contribution of the dust phase, i.e. iim sized particles, to the modification 
of the a 100 fiui layer of irradiated material is negligible, since it would only produce 
re-arrangements. For such small impactors, the cratering regime is strength-dominated. 
The volume of the crater produced Vs can be estimated as (Richardson et al 2007, Holsapple 



where Y is the effective yield strengths of the target material. We adopt a value of 
Y «i lOPa, as determined for comet Temple 1 and also in agreement with laboratory 
experiments (Richardson et al. 2007, Bar- Nun et al. 2007). Under the assumptions used 
previously, the depth of this crater for an impactor of size lO/zm is 600//m, i.e., roughly the 
penetration depth of the solar protons. 

Moreover, micro-impacts contribute less to the resurfacing because dust is carried 
away by radiation and its size distribution is much shallower than the one corresponding to 
macroscopic objects. Therefore the smallest cut-off radius that we shall assume is 1mm, 
since objects of this size are not affected by radiation. Unfortunately, there is no information 
available on the size distribution of objects of sizes ranging between 1mm and lOO's m 
orbiting in the Trojan clouds. We assume that the size distribution of the smallest observed 
Trojans (Jewitt et al. 2000) also apphes for small particles. Although this assumption may 
seem extreme, there is no creation or destruction mechanism that affects the small particles 
differently than the large objects, since none of them are affected by radiation. And the 
evidence of the cratering record on asteroids indicates that the slope of the size distribution 
of impactors is unaltered from /cm-sized bodies down to m-sized objects (O'Brien et al. 



If we take I — 1, the timescale of coUisional resurfacing is, neglecting overlapping, the 



1993): 




2006). 
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time-span in which all the surface is covered with craters, tc- Note that tc is larger than 
the time to cover the object with ejecta. Values of tq for the size range of the known Trojan 
asteroids, are shown in figure [3] for different values of the cut-off radius r^j^. For a value of 
fmin = Imm the object is completely covered by craters in about 700 to llOOyrs, while for 
Tmin = 10m, it occurs in lO'' to lO^yrs. 

Now we estimate how the area covered by optically thick ejecta, scales with the size of 
the crater. If we assume that the optically thick ejecta is close to the edge of the crater, 
and is mostly originated by material close to it, then the ejection velocity from the impact, 
Ve, scales as (Richardson et al. 2007): 



d 

Ve(x\ig -. 



Then, for the distance travelled by the ejecta, is: 

vl d 
/i oc — = -, 
9 2 



and therefore, the ratio 



h 



d/2^ 

is constant. To fix this constant we use the fact that the size of the crater produced in 
comet Temple 1 by the Deep Impact mission has a radius of approximately 30m and the 
base of the plume after the event had a radius of about 150m (Richardson et al. 2007), 
hence we adopt a value for the ratio of the areas of / = 25. In figure [3] we also plot values 
of the timescale for collisional resurfacing, r, assuming a value of / = 25. For a value of 
fmin = Imm the resurfacing occurs in about 100 to 300?/rs, while for r^in = 10m, it occurs 
in 3000 to 7000yrs. 
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3. Conclusion 

The surface properties of a Trojan asteroid are determined by the interplay of three 
different mechanisms. When an impact occurs, if it is sufficiently energetic, the inner 
materials of the body are exposed and scattered on the surface. Given that the site of 
residence and formation of the Trojans is the outer Solar System, water ice could be 
abundant, but, as shown in section [2^21 it sublimates fast from the surface, leaving a mantle 
of dust, most probably with a red spectroscopic slope. If this dust layer would remain 
unaltered for more than lO^yrs, its spectroscopic slope would turn to neutral by the action 
of solar protons. But impacts occur so frequently that the irradiation mantle is disrupted. 

The surfaces of most Trojans are red, but some neutral objects exists, therefore, it is 
suggested that the irradiation resurfacing timescale is similar, but not much shorter, than 
the collisional one. The conclusion reached here is that the timescale to fill the entire 
surface of the Trojan asteroids with craters is similar to that of the flattening of the red 
spectra by solar-proton irradiation, if the cut-off radius of the size distribution is Imm and 
its slope is —3 -for the whole size-range. If the contribution of ejecta is taken into account, 
a cut-off size of Im, with the same slope, renders a collisional resurfacing-timescale in the 
required range. 

The model presented here can be seen as an exercise to illustrate a scenario of 
collisional and radiation balance to understand the distribution of spectral slopes in the 
Trojan swarms. Naturally, refinements are needed to represent better reality. For example, 
our estimates of the collisional resurfacing timescales bear large uncertainties, mainly due 
to our lack of knowledge of the size distribution of small particles around the Lagrangian 
points of the orbit of Jupiter. On the other hand, the information regarding visual spectral 
slopes of known dynamical- family members (Fornasier et al. 2007), indicates that most 
dynamical families are remarkable uniform. With the noticeable exception of the one of 
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Eurybathes, they are mostly moderately red, suggesting a correlation with age. Eventually, 
the calculations presented here could be used to understand such a potential trend. 

Finally, a comment on how the scenario would change if the mass-loss rate of water ice 
at 120"K, based in realistic rotation rates and albedos where to be higher than estimated 
here. In this case mass loss due to sublimation would be higher and the exposed surface 
ice could have been lost from the Trojans in their current orbit. A different orbit is not 
required. 
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Fig. 1. — Histogram of the slopes of visual reflection spectra normalized at 6OOOA for the 
Trojan asteroids. The sources of this data can be found in Melita et al. 2008. 
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Fig. 2. — Slopes of visual reflection spectra normalized at 6000A vs. D(0.04), the physical 
diameter, assuming an albedo of 0.04, for the Trojan asteroids. Data taken form Melita et 
al. 2008. 
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Fig. 3. — Left: Timescale in which the whole body is covered with craters, tc, as a function 
of the radius of the target, for the size range of the known Trojan asteroids. Right: Timescale 
of coUisional resurfacing, assuming a value of i = 25, r, as a function of the radius of the 
target, for the size range of the known Trojan asteroids. See text for details. 
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